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ABSTRACT 

The distance to the black hole binary GX 339-4 remains a topic of debate. We examine high-resolution 
optical spectra of the Na D lines resolving the velocity structure along the line of sight. We find this 
to be complex, with at least nine components, mostly blue-shifted, spanning a velocity range of nearly 
200 km s -1 . The presence of components with a large blue-shift rules out a nearby location and 
requires that the binary be located at or beyond the tangent point, implying a lower limit to the 
distance of ~ 6kpc. The presence of a significant red-shifted component at -fSOkms^ 1 is even more 
intriguing as GX 339-4 also has a slightly positive systemic velocity, suggesting that the source, and 
this cloud, could be on the far side of the Galaxy, where the radial velocities due to Galactic rotation 
become positive again. If this is the case, we require a distance of ~ 15kpc. This is less secure 
than the 6kpc lower limit however. We discuss the implications of these possible distances for the 
outburst and quiescent luminosities, and the nature of the companion star, and argue that a large 
distance is consistent with these characteristics. In particular, it would explain the non-detection of 
the companion star during the faintest states. 
Subject headings: stars: individual: GX 339-4 



1. INTRODUCTION 

The distance of X-ray binaries is challenging to deter- 
mine precisely. They are usually too distant for current 
parallax measurements, and there are no really reliable 
standard candles, other than type-I X-ray bursts show- 
ing radius expansion, and the se only occur in a su bset 
of neutron star systems (see iKuulkers et alJ 120031 and 
references therein). For black holes, such as GX339- 
4, the problem is much less tractable. Fortunately 
the majority of such systems undergo periods of deep 
quiescence when it is possible to detect the compan- 
ion star, and determine the binary parameters. Hence 
both the radius and temperature (i.e. surface bright- 
ness) can be estimated and th e distance derived from the 
appa rent brightness (see e.g. lOelino. Harrison. & Orosd 
l20f)l . For GX 339-4 this is not possible, as even in a 
very low luminosity state, no spectral features of the 
companion star were seen in VLT /FORS observations 
(Shah baz. Fender, fe Charles! 120011) . This is frustrating 
for distance and parameter estimation, and very puzzling 
in its own right. One explanation for the non-detection of 
the companion is that the distance has been significantly 
underestimated. 

Current distan ce estimates are very uncertain; 

iBuxton fc Vennesl l)2003l) has provided a convenient sum- 
mary. Early estimates were rather low, mostly falling in 
the 1-4 kpc range. More recently, larger estimates have 
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been proposed. The non-detection of the companion star 
( Shah baz~Fender. fc Charlesl2001j) was used as evidence 
for a dis tance >, 5.6 kpc. A more novel approach was 
taken bv iMaccarond l)2003|) who has shown that in most 
black hole binaries, transitions between the high/soft and 
low/hard states occur at approximately a fixed fraction of 
the Eddington limit. It was suggested that the transition 
luminosity can then be used as a mass-dependent stan- 
dard candle. Assuming a minimum black hole mass for 
GX 33 9-4 of > 5.8 M Q ijHvnes et a.l.|l200l . IMaccaronel 
(2003) estimate a distance > 6 kpc, allowing for scatter 
about the relation. 

The actual distance of GX 339-4 thus remains very 
uncertain, but could be rather large. This introduces 
considerable uncertainty into interpreting its quiescent 
and outburst luminosities, and the nature of the com- 
panion. The distance also has important implications 
for interpretation of the recently discovered resolved jets; 
I G alio et al.l l)2004fl estimated that the motion of the jet 
head requires an apparent velocity of 0.9 c even at 4 kpc. 
A significantly larger distance would therefore imply ap- 
parent superluminal motion. We have therefore analyzed 
interstellar features in high resolution spectra described 
in Section[2]to elucidate the properties of the line-of-sight 
in more detail. In Section [21 we perform a detailed study 
of the properties of line-of-sight Na 1 and Ca 11 absorption, 
resolving most individual clouds in Nai and hence de- 
termining cloud-by-cloud velocities, velocity dispersions, 
and Nai column densities. We also have some informa- 
tion about the relative abundances of Can to Nai as a 
function of velocity. We discuss the implications for the 
distance of GX 339-4 in SectionfU and for the reddening 
in Section [SJ Finally in Section we consider how the 
range of distances proposed affect estimates of the out- 
burst and quiescent X-ray luminosities and the nature of 
the companion star, and in Section [7] we summarize our 
conclusions. 
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2. OBSERVATIONS 

High- resolution observations of GX 339-4 were ob- 
tained with the UV- Visual Echelle Spectrograph (UVES) 
on UT2 (Kueyen) at the VLT on 2002 August 9-15. 
These spanned a wavelength range of 4100-6200 A and 
included a total exposure time of 12 hr yielding a superb 
quality average spectrum (signal-to-noise ratio near the 
Na D lines greater than 100 without binning). A 1" slit 
was used giving a nominal resolution near the Na D lines 
of ~ 44 000, corresponding to a kinematic resolution of 
6.8kms~ 1 . The starlight will be concentrated to the cen- 
ter of the slit, however, yielding a somewhat higher reso- 
lution. From the spatial profiles near Na D we estimate 
a mean seeing of about 0.5" , implying that the resolution 
could be as good as 4kms _1 with a 1" slit. In practice, 
this will not be realized and we expect a value between 
these extremes. We therefore leave the instrumental res- 
olution as a free parameter in subsequent fits. For the 
rest of this work we will focus only upon the Na Dl and 
D2 lines (5895.92 A and 5889.95 A respectively). Pipeline 
optimal extractions were supplied and were of good qual- 
ity in this region; a manual reduction of selected spectra 
yielded negligible improvement. 

We supplement these data with lower-resolution ob- 
servations of the Call H and K lines (3968.47 A and 
3933.66 A respectively). These were obtained with the 
RGO spectrograph of the Anglo- Australian Telescope 
over 2002 June 6-11, with the R1200B grating yield- 
ing a wavelength coverage of 3500-5250A. The images 
were de-biased and flat-fielded, and the spectra subse- 
quently ext racted using conventional optimal extraction 
techniques l|Hornei il986'l. Wavelength calibrations were 
interpolated between CuAr comparison lamp images, ob- 
tained every 20-30 mins. The kinematic resolution is es- 
timated to be 74kms _1 . 

3. LINE-OF-SIGHT VELOCITY STRUCTURE 
3.1. Nal absorption 

We show line profiles of the two Na D lines in Fig.Q]rel- 
ative to the local standard of rest (LSR) . They are clearly 
complex with a large velocity dispersion, suggesting a rel- 
atively large distance. This immediately rules out sug- 
gestions that GX 339-4 is a local object at ~ 1 .3 kpc 
l|Mauche fc Gorensteinl 119861: [Predehl et all 119911): the 
same conclusion was drawn bv lZdziarski et alJ l)1998fl for 
different reasons. 

We can obtain an approximate estimate of the distri- 
bution of Na I with velocity by converting the line pro- 
file s into apparent optical depths (Fi g. HI), as discussed 
by iSembach. Danks. fc Savaga l) 19931) . This is not the 
true optical depth, because it is modified by instrumental 
broadening. It is also reduced (relative to the Nal col- 
umn density) for strongly saturated components. This 
is seen in differences in the optical depths inferred from 
the Dl and D2 lines. We have scaled- up the Dl op- 
tical depths in Fig. [21 by a factor of two (the ratio of 
oscillator strengths) so that the two lines can be directly 
compared. Where the two optical depth profiles over- 
lay (i.e. the Dl line has exactly half the optical depth 
of the D2 line) the apparent optical depths are simply 
blurred versions of the true optical depth, and the rela- 
tive contributions of different components are preserved. 
Where unresolved saturated structure is present, the D2 
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Fig. 1. — Line profiles of the Na Dl (dashed) and Na D2 (solid) 
lines. Spectra have been normalized, but no offset is applied be- 
tween the lines. Labels denote distinct components. Some of these 
(e.g. D) appear to be further unresolved multiples. Tick marks 
indicate the components fitted in Table HI 
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Fig. 2. — Apparent optical depths of the Na D lines. The D2 
optical depths have been doubled so that they are identical to Dl 
in the optically thin limit. Where the lines disagree, saturation is 
important and the optical depth may be underestimated. 



line becomes preferentially more saturated, the Dl line 
becomes proportionately stronger, and the ratio of ap- 
parent optical depths approaches unity. In this case the 
total absorption of a component is underestimated, but 
the Dl optical depths reflect the column densities more 
accurately. We can see that most of the absorption col- 
umn appears to be at low velocities, |Vlsr| 20kms _1 , 
especially when one considers that these are the most 
saturated components, and hence the most severely un- 
derestimated. There is also a significant contribution at 
higher velocities, however. 

We can perform a more quantitative analysis by mod- 
eling the line profiles, and explicitly accounting for un- 
resol yed saturation. We follow standard t echniques 
(e.g. lSpitzerill978HSembach. Danks fc Savagdll993T) . as 
also applied bv iDubus et alJ l )2001[) to Call lines in 
XTE Jl 118+480. We assume that the line-of-sight ab- 
sorption consists of a number of discrete clouds, each 
of which has mean velocity v, a Gaussian distribution 
of velocities with width 6, and Nal column density n. 
Given this model we can simultaneously predict both Dl 
and D2 intrinsic line profiles. We then convolve these 
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with the instrumental resolution and compare with the 
data. As discussed above, the exact instrumental resolu- 
tion is uncertain and was fitted as a free parameter. The 
best fit was obtained for 6.0 kms , which falls within 
the range expected given the slit width and seeing (Sec- 
tion[3J). Note that in most cases these parameters are rel- 
atively well constrained, even though the intrinsic widths 
are not fully resolved. This is because the differential sat- 
uration of the Dl and D2 lines makes the ratio of the two 
line strengths sensitive to the intrinsic widths. Hence we 
can recover a reasonable estimate of the column density 
even for somewhat saturated lines. The largest uncer- 
tainty exists where significant saturation is present, and 
then we sometimes only constrain a lower limit upon the 
column density. Fitting was performed iterative ly us- 
ing t he Levenberg-Marquardt algorithm (see e.g. IPressI 
1988). We began by fitting the components labeled in 
Fig. ^ Additional components were added when it was 
clear that they were needed; for example D clearly in- 
cludes an additional low velocity component and a weak 
component is also needed between G and H. Other defi- 
ciencies were apparent only by comparison with the pre- 
liminary fits and lead us to add additional components 
near C and F as well. It may be that such weak compo- 
nents from small clouds are actually present throughout 
the profile, but are only discernible in the gaps between 
the larger components (Redfield priv. comm.) 

The derived parameters for the components are listed 
in Table HI This model reproduced the data extremely 
well, with no significant residuals, so we have not shown 
the fit. Velocities are derived from both Dl and D2 lines 
simultaneously, and indicate a range of mean cloud ve- 
locities of -143 to +32 kms -1 . We will compare this 
to Galactic rotation curves in Section Component 
widths a re l.S-T.Skms -1 . These values are typical of 
the ISM: ISemha,ch Danksl flflfll found average values 
for the Na D lines of 4.0 kms -1 . If the lowest widths 
are purely thermal, they indicate cloud temperatures of 
> 3 000K, but turbulent broadening could be present as 
well. It may be that some of the larger widths, such as 
for component B, indicate additional unresolved compo- 
nents rather than intrinsically high widths. One of these 
mechanisms is needed for the highest widths, as these 
are too high for thermal broadening alone. We derive 
a total Nai column density of at least 3.5 x 10 13 cm -2 , 
of which components B and C (assumed to be relatively 
local) contribute at least 63%. 

In most cases the sensitivity to unresolved blends is 
weak as the components are not strongly saturated and 
the total column density is approximately conserved. For 
component B, however, the total column density is sen- 
sitive to assumed substructure. To illustrate this, we 
repeated the fitting with two sub-components; the re- 
sults are shown in Table This allows the main com- 
ponent to be narrower, and hence saturate more eas- 
ily. A much larger column density is then required to 
achieve the same equivalent width of absorption. Con- 
sequently, the single component fit should be considered 
the lower limit on the column density in the BC complex 
(2.2 x 10 13 cm -2 ). This implies that the total column 
density estimated above, and the fraction of total column 
contained within the BC complex, are also lower limits. 
Similar difficulties exist for components CI and HI, and 
large column densities cannot be formally rejected in ei- 



Table 1 

Results of fits to the Na D lines. 



Component 



v a b h 
(kms" 1 ) (kms" 1 ) (10 1 



Single component B 
32.4;' ' 
25.4" 



Al 
A2 
B 
CI 
C2 
Dl 
D2 
E 
Fl 
P2 
G 
HI 
H2 
J 



(1.3 
1-1.2 
-0.4 

3.0 ±0.1 



-13.7 
-23.0 



+0.2 
-0.3 
+ 1.8 



-35.4 ±0.3 2.2 



1 5+ ' 2 
9 q+2.3 
z - z -0.5 

7.3 ±0.2 

7+0.5 
2.7 
+ 1.9 
-2.4 
+0.5 
-0.3 



2.7] 
7.4 



0.80±p 

o.45tg;« 

12 7+°' 7 

iz -'-0.5 
».3_ 2 6 

1.17 ±0.35 



0.94 



-46.0 ±0.1 4.0 ±0.2 3.05 
-65.8 ±0.1 5.1 ±0.2 2.04 ±0.05 



+0.10 
0.09 
+ 0.14 
-0.12 



-82.2 



+0.3 



2-l_ 
6.6^ 



+0.6 



-0.4 

-90.lti;° 

-106.8 ±0.2 3.7±o| 

+0'.6 
-0.4 



0.53 ±0.08 
45+ 07 

u '^°-0.08 

1.10 ±0.04 
0.44 042 



-117.5 ±0.4 1.2_ 4 v.«_ . 05 
-128.6 ±0.1 4.4 ±0.3 2.32±°°l 

7+0.4 q 75 +0.'06 



-143.0 ±0.2 2.7 H 



0.3 



0.04 



Split component B 
Bl 12 -6lo.7 
B2 2.3 ±0.1 

Cl c -13.3 ±0.2 5.5 



+0.7 



C2 



-26.4 



+0.6 

-0.5 



^•■'■-o.s 

3 0+ - 4 

+0.4 

0.5 

3 1+ 1 ' 4 
g^-iq 



1.06 



+0.21 
-0.13 
1060 



250l 150 

-,+0.16 
-0.15 
53+° 12 



5.20"* 



a Mean velocity. 
b Gaussian width. 
c Nal column density. 

d Revised parameters for B and C. Other components are unaf- 
fected. 

c There is no formal lower limit on the width of component CI, or 
a corresponding upper limit on the column density. 



ther case. This reflects the significant saturation also 
present in these components (Fig. |2J • 

3.2. Call absorption 

The resolution of the AAT observations of the Can 
lines is not sufficient to perform a similar component-by- 
component analysis. We can, however, use the model de- 
rived from the Na D lines to predict the lower-resolution 
profiles and perform a crude comparison. To do this we 
initially assume a ratio of Call to Nai column densities 
of 0.6, appropriate for in ner Galactic center sight-lines 
(Scmbach & Danksl fl994|) . and an instrumental resolu- 
tion of 74 km s -1 . The intrinsic width of each component 
is scaled down by 0.76, which is appropriate if it is dom- 
inated by thermal broadening (i.e. scaling by the square 
root of the atomic weight). The fit to the profile of the 
K line is shown in Fig. |2| the H line is not used because 
it is blended with Balmer emission, but the line profile 
otherwise appears to be similar. The abundance ratio 
assumed clearly underestimates the Can absorption at 
higher velocities. To improve the fit, we recalculated as- 
suming a small velocity offset (5 kms -1 ; a small fraction 
of a resolution element) for the AAT spectra, a slightly 
lower ratio of Call to Nai of 0.53 for most components, 
but a higher ratio of 4.5 for components G-J. If we al- 
ternatively keep the intrinsic width, 6, fixed (appropri- 
ate for turbulent broadening) then we require abundance 
ratios of 0.42 and 2.5 for low and high velocity compo- 
nents respectively. Either model provides a much bet- 
ter fit, although clearly still underestimates the wings of 
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Fig. 3. — Line profile of the Can K line. Histograms show 
the data, at much lower resolution than for the Na D lines. The 
dashed line shows a 'blind' prediction of the line profile based on 
the Na D fits, assuming no velocity correction is needed, that the 
intrinsic component widths are thermal, and that the ratio of Ca II 
to Nal is 0.6, appropriat e for inner Galactic center lines of sight 
ISembach fc Danks 199|). The solid line is a better fit including a 
small velocity offset and allowing the abundance of the G— J com- 
plex to differ from the rest of the components. See the text for 
more details. 



the Can profiles. This is a common feature; additional 
Can absorption is often seen at high velocities (beyond 
those allowed by Gala ctic rotation) where Nal absorp- 
tion is weak or absent t'Scmbac h fc Dankslll994|) . It has 
been attributed to a significantly enhanced abundance of 
gas-phase calcium (by a factor of 50 or more relative to 
sodium) in the high velocity components of the ISM via 
collisional destruction of dust grains. It is also observed 
that the m ean calcium abundance is enhanced along halo 
sight lines l)Sembach fc Dankslll994(L although this may 
simply reflect the preponderance of high velocity mate- 
rial at high latitudes. Other than such diagnostics, the 
Nal lines are more useful, as they are more closely related 
to Galactic rotation, and correlate better with E(B — V). 

4. THE DISTANCE TO GX 339-4 

We have mapped out the properties of gas clouds along 
the line-of-sight to GX 339-4 in velocity space. We now 
attempt to associate these with distance, assuming the 
velocity is due to Galactic rotation, and hence derive a 
lower-limit on the distance to GX 339-4. Qualitatively 
we expect to see low-velocity material from the Solar 
neighborhood, and increasingly negative velocities as the 
tangent point is approached. Beyond this, velocities be- 
come less negative again, and become positive at dis- 
tances outside the Solar circle. The preponderance of 
absorption components with negative velocities is thus 
as expected from Galactic rotation along this sight line. 
To be more quantitative, we show the expected velocity 
as a function of distance in Fig.lH We assumed the rota- 
tion curve prescription of Clemens ( 1985) , with the mod - 
ified polynomial coefficients of Nakan ishi fc Sofud |2003). 
These were derived for a Solar Galactocentric radius of 
i?o = 8.0 kpc and velocity of 6o = 217 km s _1 , consis- 
tent with the most recent geometric determination of 
the distance to the Galact ic center (Rp = 7.94±0 .42kpc, 
9 = 220.7±12.7kms- 1 : lElsenhauer et all2003lK Fig.|l 
also indicates the velocity and relative column density of 
each fitted absorption component for easy comparison 
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Fig. 4. — Expected relation between radial velocity and distance 
along the line of sight to GX 3 39-4. The c urve uses the rotation 
curve adopted by Nakanishi &: Sofud 120031) . Each horizontal bar 
represents one of the fitted components of Na D absorption. The 
vertical position indicates the fitted velocity and the length of the 
bar is proportional to the column density. The horizontal positions 
are arbitr ary. The hatch e d reg ion indicates the distance range 
favored by Zdziarski ct al. 11998). 



with the predicted velocities. 

Components B and C in GX 339-4 likely correspond 
to local material, and contribute > 63 % of the to- 
tal Nal column density. The velocities correspond to 
a distance < 2 kpc, but the material could be within 
1 kpc allowing for a lOkms -1 dispersion in mean ve- 
locities. For a Galactic latitude of —4.5°, this mate- 
rial will thus be < 0.16 kpc below the plane and well 
within the Na I layer (which has an estimated scale height 
0.43 + g;po kpc: ISembach fc DanksllT99l . Examination of 
field stars near GX 339-4, at distances up to 2 kpc, sug- 
gests that local material within this distance accounts for 
~ 2/3 of the reddenin g to GX 339-4, consistent with this 
HZdziarski et alJlTflfl^. 

If the distance to GX 339-4 were ~ 4 kpc, as pro- 
posed by iZdziarski et alJ (^998) and iBuxton fc Vennesl 
(2003), then the maximum expected negative velocity 
would be ~ — 60kms~ 1 . Components D and E could 
be co-rotating, but the G-J complex would require pe- 
culiar velocities of 40-80 km s -1 . The latter accounts 
for < 13% of the total column density, and as we 
have already noted is relatively overabundant in Can. 
There are two possible interpretations for this high ve- 
locity absorption. We consider the most likely to be 
that it is co-rotating with the Galactic rotation curve. 
In support of this, the similarity of the velocities of 
the G-J complex to the predicted tangent point veloc- 
ity of approximately — 120kms is suggestive. The 
range of velocities is also in very good agreement with 
CO maps, which exhibit a vertically integrated range 
of velocities to — 140kms~ 1 for this Galactic longitude 
ijDame. Hartmann. fc Thaddeusll2001[) . The inferred al- 
titude of the tangent point material, is ~ 0.55 kpc. This 
is not a serious problem given a scale height of 0.43 kpc. 
The alternative explanation is that the G-J complex 
(and component F) arises from closer material with large 
peculiar velocities; the same materia l as th e 'forbid- 
den velocity' gas of ISembach fc Danksl |[l994). The en- 
hanced Can abundance would be typical of such gas. 
We believe that the presence of multiple components 
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contributing a significant fraction of the Nai column 
at such high velocities is very unlikely, however. Such 
features are not see n in Nai in any of the similar line s 
of si ght collected by ISembach. Danks. fe Savage) l)1993j) . 
and ISembach fc DanksT " ifT994l) estimated that even the 
velocity dispersion of the high velocity Can absorption 
components was only 21.3kms _1 . If GX339-4 lies at 
4kpc it would therefore have a rather unusual line-of- 
sight. Neither are such extreme velociti es seen in any of 
the sp ectra of local stars presented by iWeltv fc Hobbsl 
(2001). Consequently, it seems more likely that most 
of the absorption velocities do indicate Galactic rotation 
and that GX 339-4 lies at or beyond the tangent point. 
The implied distance is then > 6 kpc, allowing for some 
uncertainty in rotation curves and the intrinsic velocity 
dispersion of the absorbing clouds. 

Beyond this, the situation is less clear. This is be- 
cause the velocities increase again, and so there is am- 
biguity between close and distant material. Compo- 
nents D-F could arise on either side of the tangent 
point, and it is hard to make a conclusive statement 
about how far beyond this GX 339-4 could lie. Two 
lines of argument suggest a significantly larger distance 
is possible, but neither are conclusive. One is the sys- 
temic velocity of GX 339-4 itself. The best estimate of 
this comes from the emission line r adial velocity ana l- 
ysis of the same dataset used here ijHvnes et all 12003). 
with Hen emission line wings suggesting 7 = +12.3 ± 
l.Okms - 1 , and Bowen lines consistent with this for the 
preferred orbital period. This would of course be ex- 
pected for a local object, but the Na D velocity dis- 
persion rules this out, even if earlier arguments do not. 
If GX 339-4 lies at the tangent point, however, then 
the peculiar radial velocity would be ~ +140 km s -1 . 
This is unusually large for a black hole binary, as 
iWhite fc van Paradiisl l) 19961) estimated that their ra- 
dial velocity dispersion should be 40kms~ 1 , based on 
their distances from the plane and consistent with 
measurements. Larger values are not impossible (e.g. 
GRO J 1655-40 h as an estimated peculiar radial veloc- 
ity o f 114kms _1 ; iBrandt. P odsiadlowski. & Sigurdsson 
1995), but are thought to arise when black holes are 
formed by fallback onto a neu tron star which received 
a kick l|Frver fc Kal ogera 20(H1). A large velocity is only 
expected for a low-mass black hole, however, which is in- 
consistent with the high mass function of GX 339-4 and 
lack of evidence for a high inclination. Further increasing 
the distance is therefore appealing as it would reduce this 
required peculiar velocity. A location close to the Solar 
circle on the far side of the Galaxy (> 15 kpc) would 
require negligible peculiar velocity, consistent with most 
Galactic black hole binaries. The other argument comes 
from the one absorption feature we have not yet ac- 
counted for: component A. This has two sub-components 
with positive velocity, 25-35 km s _1 . If this arises locally 
then it arises from material with a forbidden velocity 
> 25kms _1 . Given that there are only two weak com- 
ponents (or one with a complex profile), and that the 
total contribution to the column density is < 4%, this is 
not impossible. This material could also arise on the far 
side of the Galaxy, and would then have a radial velocity 
close to that of the source. One could argue that since 
GX 339-4 has a Galactic latitude of —4.5°, at this dis- 
tance, the line-of-sight should be more than 1 kpc below 



the plane. However the Galaxy does exhibit warps, and 
the outer d isk warps downward towa rd the line-of-sight of 
GX 339-4 (jNakanishi fc Sofud 12003]) . It is possible that 
GX 339-4, and the red-shifted absorption, are associated 
with this warped material. 

We thus conclude that a minimum distance of 6 kpc is 
required, based on the presence of intervening material 
with large blue-shifts, consistent with expected tangent 
point velocities. The presence of red-shifted absorption, 
together with the low systemic velocity of GX 339-4 it- 
self, suggest that a much larger distance > 15 kpc might 
be correct, although a large peculiar velocity for GX339- 
4, together with forbidden velocity gas could obviate this. 
GX 339-4 is thu s further away than most previous es- 
timates (but see iShahbaz. Fender, fc Charles! 120011 and 
lMaccaronell2003D . and possibly a lot further. 

5. THE REDDENING OF GX 339-4 

Our high quality spectra should also allow us to obtain 
more robust reddening constraints from the Na D lines; 
we do not attempt this wi th the Call lines as thes e do 
not trace the dust as well l)Sembach fc Danks 1994). 

The integrated Na D equivalent width (EW) has been 
estimated by a number of authors. The most recent val- 
ues are 3.7 A llSoria. Wu. fc Johnstor]IT999T) and 4.2 ± 
1.0 A ijBuxton fc Vennesl l2003|) . We estimate a total 
EW of 3.57 ±0.02 A from our resolved Na D lines, with 
1.48 ± 0.01 A and 2.10 ± 0.01 A in the Na Dl and D2 
lines respectively, consistent with other estimates. The 
uncertainties quoted are conservative estimates based 
on the range of possible continuum levels; formal sta- 
tistical errors are smaller than this. The interpreta- 
tion of the equivalent width is obviously complicated 
by saturation of some components, but neglecting this 
we recover a lower-limit of E(B — V) > 0.6 from the 
Na D2 line assuming the asymptotic linear calibration of 
iMunari fc Zwitterl l)1997D . We can perform an analogous 
calculation using the less saturated Na Dl line, since in 
the optically thin limit the D2/D1 ratio has a value of 2. 
This yields a lower limit of E(B -V)> 0.85. 

An alternative approach is to circumvent the satura- 
tion effects by using the column density derived from 
the line profile fits. As with all such calibrations, it is 
still sensitive to assumptions about gas-to-dust scaling, 
ratios . We use the values derived bv ISembach fc D anks 
(1994), which include scalings for a number of types of 
sight line. An Nai column of at least 3.5 x 10 13 cm -2 cor- 
responds to A H = 3.1 x 10 21 cm- 2 and E(B - V) = 0.50 
using ratios appropriate to inner Galactic center sight 
lines. These values change to Ah = 3.3 x 10 21 cm -2 
and E(B — V) = 0.69 using Galactic center ratios, and 
A H = 4.4 x 10 21 cm" 2 and E(B - V) = 0.59 for av- 
eraged ratios. The uncertainty is larger than implied 
by the scatter between these estimates, however, since 
the variance in ratios within each of these samples is 
large. While some of the variance may simply reflect 
measurement errors, it is likely that the ratios are sub- 
ject to significant cloud-to-cloud scatter. As discussed 
earlier, these are also lower limits, since the local compo- 
nents of Na D are saturated and the column density used 
here may have been significantly underestimated. This is 
likely why both the Ah and E(B — V) val ues we estimate 
are lower than other measurements (e.g . [Zdziarsk i et all 
ll998tlCorongiu et aUEffll IMiller et alll2004lK 
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Consequently, even with the detailed resolution pre- 
sented here, we cannot determine a reliable E(B — V) 
from the Na D lines; this method is fundamentally rather 
crude, and limited by saturation. All we can say with 
confidence is that the reddening is greater than 0.5, and 
probably greater than about 0.8. This i s comparable to 
previ ous estimates i ncluding 1.2 ± 0.1 llZdz iarski etahl 
fl998^ and 1.1 ± 0.2 {Buxton fc Vennesll2003|) . We sus- 
pect that these quoted uncertainties are underestimates, 
however, as similar difficulties apply as with our own 
measurements; indeed the latter estimate is also based 
on Na D. Upper limits are probl emat ic, and neither the 
metho ds of lZdziarski et al.1 (0,998) , nor lBuxton fc Vennesl 
( 2003), nor our discussion above adequately address this. 
Neutral hydrogen column densities do allow a crude up- 
per limit, however. A range of values are derived (see 
references above), but all imply Ah < 7 x 10 21 cm~ 2 . 
Similarly, a range of gas-t o -dust scalings exist, but that of 
IBohlin. Savage. fc Drakel ©78), with an assumed 30% 
uncertainty, encompasses all of them, and also accounts 
for the variance between lines of sight. Using this scaling 
and the upper limit on Ah implies E{B — V) < 1.6. 
The actual value could well be much less than this. 
An inde pendent check is provided by t he Ga lactic dust 
maps of Schlegel, Fink bemer. fc Davisl {19981) which im- 
ply E(B — V) ~ 0.9, although their use for Galactic lat- 
itudes less than 5° is not recommended by the authors. 

6. DISCUSSION 

6.1. The outburst luminosity 

The luminosity of the transition between hi gh/soft and 
low /h ard states in black holes is of interest (fMaccarond 
12003]) and it may even b e an approx i mate standard can- 
dle. As already noted, IMaccaronel lj2003j) estimated a 
minimum distance of 6 kpc assuming the minimum black 
hole mass is given by the mass function of 5.8 ± 0.5 M 
<|Hvnes et alJ l20031. 'The true mass is likely to be signif- 
icantly larger than this, however, as the mass function is 
only approached for high inclinations and there is no evi- 
dence that GX 339-4 has a high inclination; in deed it has 
been argued that the inclination is very low l)Wu et al.l 
12001(1 . For actual black holes masses of 10, 15, or 20 M , 
distances of about 10, 12, or 14 kpc respectively woul d 
be expected based on the method of IMaccaronel {2003). 
These are very reasonable if GX 339-4 lies somewhat be- 
yond the tangent point as we have argued. If it is really 
beyond 15 kpc, however, then either the black hole mass 
is relatively high, or the transition luminosity is rather 
h igh (but not imp l ausib ly so). 

Zdziars ki et al.l (|1998j) have instead tried to use the 
maximum observed luminosity as a constraint on the 
distance. We can reverse the argument, however and 
examine the maximum luminosity for different possible 
distances. At 6 kpc, their quoted maximum observed 
luminosity corresponds to just 10% of the Eddington 
limit of a 10 M Q black hole. This seems rather low, 
and a higher black hole mass would further exacerbate 
this. In contrast, for a distance of 15 kpc and black 
hole masses of 10-20 M Q , the Eddington ratio is 0.3- 
0.6, which seems much more in line with other Galactic 
black hole X-ray binaries. More recent observations pro- 
vide more stringent constraints, yielding an 0. 1-200 keV 
flux of 5.4 x 10~ 8 ergcm~ 2 s _1 (extrapolated from a 3- 
100 keV fit; Homan et al. in preparation). For a 10 M 



black hole at distances of 6 and 15 kpc respectively this 
corresponds to 20% or 110% of the Eddington limit. 
The latter could readily be reconciled with a black hole 
somewhat more massive than 10 M Q . Thus the luminos- 
ity is consistent with a relatively large distance given the 
implied high mass for the black hole. 

6.2. The 'quiescent' luminosity 

The X-ray luminosity in the lowest observed state 
( Galk xF^ender. fc Corbell |2003) corresponds to 2 x 
10 33 ergs _1 or 1 x 10 34 ergs _1 for distances of 6 or 
15 kpc respectively. Most quiescent black holes have 
luminosities of between 2 x 10 30 and 6 x 10 31 ergs -1 
llGarcia e t alJl200 lt ISutana et al.ll2002t iHa meury et alJ 
l2003HTomsick et al.ll2003UMcClintock et alll2003l). The 
excep tion is V404 Cyg (5 x 10 33 ergs~ 1 : iGarcia et alJ 
I2001I) . This is a long period object (6.5 day) and may 
be something of an anomaly among quiescent black 
holes; indeed it exhibits dramatic X -ray variability 
{Wagner et alJ 11994 iKong et all 120021: Hynes et al. in 
preparation) which suggests that it may not be in a fully 
quiescent state. GX 339-4 appears to have a compara- 
ble or larger luminosity at minimum light, so its lowest 
observed state may also not correspond to the typical 
quiescent states of transient black hole binaries. 

6.3. The nature of the companion star 

A relatively active system is appealing in the light of 
the non- detection of spectral features from the compan- 
ion star (Shah baz. Fender, fc Charlesll2001jL Two expla- 
nations exist for this: either the accretion flow is so lumi- 
nous that it masks the companion star, or the companion 
is an early- type star with no spectral features in th e re- 
gion studied bv lShahbaz. Fender, fc Charlesl {2001D . As 
noted above, a relatively active accretion flow already 
favors a larger distance, and even 15 kpc gives a plau- 
sible luminosity. If we require an early-type companion 
then its luminosity will be larger than for the typical K 
or M companions of Galactic black holes, and so a large 
distance might also required to avoid it being visible. 

With no detection of the photospheric lines of the com- 
panion, its nature is poorly constrained. The companion 
star s hould be somewha t evolved, as a 1.7 day orbital pe- 
riod {Hvnes et al J 120031) implies a mean density for the 
companion of only 0. 06 g cm" 3 . This is a factor of t en be- 
low that obtained by Shahb az. Fender, fc Charlesl {200 If) 
who assumed a shorter orbital period. Such a low den- 
sity is reached on the main sequence only for B0 stars; 
it is unlikely that the companion is such a star, however, 
as it would then be expected to dominate the total light 
at low luminosities, and also significantly affect the mo- 
tion of the black hole. Hence the companion of GX339- 
4 has a lower density than a main-sequence star of the 
same temperature. It is probably a sub-giant, as late- 
type giants have even lower densities than we derive for 
GX 339-4. 

For relatively early types, the companion mass is also 
significant and should res ult in measurable o rbital mo- 
tion of the accretion disk. iHvnes et alJ {2003D found no 
such evidence (only a weak modulation was present, with 
the wrong phase-dependence) and suggested a mass ra- 
tio of q < 0.08. This argues for a companion mass 
of 1.6 M© or less unless the black hole mass is signifi- 
cantly larger than in other objects (> 20 M ) and hence 
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a spectral type of F or later unless the companion is 
hotter than a main sequence star of the same mass. 
In support of these arguments, GROJ1655-40 has a 
comparable orbital period (2.6 d), and an F companion 
l|Orosz fc Bailvnll997tlShahbaz et alll999ft . Its mass ra- 
tio is 0.42 ± 0.03 l|Shahbazll2003D . however, t he disk lines 
move with an amplitude of 76.2 ±7.5 km s _1 ijSoria et alJ 
I1998J) . and the companion totally dominates the optical 
light in quiescence and remains visible even during out- 
burst. These properties are very different to GX 339-4, 
which evidently has a much less influential companion. 
G or K subgiants therefore seem more likely than earlier 
types. 

At our proposed distances there is little difficulty in 
accommodating such companions. We consider a num- 
ber of examples. For each we adopt the main sequence 
mass corresponding to the spectral type, but adjust the 
radius to fill the Roche lobe. The size of the lobe is 
constrained by the 1.7 day orbital period, an assumed 
black hole mass between 5.8 and 20 Mq, and a mass ra- 
tio less than 0.1. The reddening is assumed to be be- 
tween 0.8 and 1.6, as discussed above. The upper limit 
is more pertinent, as it makes it easier to hide luminous 
companions. An F0 star (adopted mass 1.6 M ) could 
only be accommodated if the black hole mass were near 
20 M©; otherwise we would expect detectable motion of 
the accretion disk. At 6 kpc such a star should have an R 
magnitude of 17.0-19.5, dependent on reddening. This is 
brighter than the observed faintest magnitude (R = 20.1; 
Shah baz" Fender, fc Gharle"s!l200lD . so is unlikely unless 
the companion is significantly less massive (and hence 
has a smaller Roche lobe). At 15 kpc, however, this could 
be as faint at R — 21.5, so could be masked by the ac- 
cretion light. For a GO star (assuming M — 1.1 Mq), 
we expect a range of 18.1-20.6 at 6kpc, or 20.1-22.6 at 
15 kpc, and either distance is acceptable, provided the 
reddening is relatively large. A K0 star (M = 0.8 Mq) 
yields R = 19.1 - 22.7 at 6 kpc and 21.1-24.7 at 15 kpc 
and poses no difficulties. 

These estimates are rather crude and dependent on 
some assumptions, but are probably still useful as ap- 
proximations. To summarize, the companion to GX339- 
4 is likely to be a sub-giant, of spectral type G or later 
if close, ~ 6 kpc, or possibly F if the distance is much 
larger. 

7. CONCLUSIONS 

We have analyzed Nai and Can absorption along the 
line of sight to GX 339-4. We find a rich absorption spec- 
trum with a large velocity dispersion. Velocities are pre- 
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dominantly blue-shifted and span the range of velocities 
expected for Galactic rotation along this line-of-sight. 
The presence of multiple features with a significant col- 
umn density near the tangent point velocity suggests that 
GX 339-4 is either at or beyond the tangent point, i.e. 
d > 6 kpc. The presence of absorption at red-shifted ve- 
locities, together with the systemic velocity of GX 339-4 
itself, suggest the possibility that the object could be at 
an even greater distance, on the far-side of the Galaxy, 
d > 15 kpc. Neither of these pieces of evidence are as 
secure as that for a location beyond the tangent point, 
however. 

Our analysis also indicates a total Na I column density 
along the line of sight of at least 3.5 x 10 13 cm~ 2 , and this 
suggests a reddening of at least E(B — V) = 0.5, depend- 
ing on assumed ratios of Nai to dust densities. This is 
consistent with more straightforward calibrations of the 
equivalent widths, which imply E(B — V) > 0.85, but it 
could be considerably larger dependent on the degree of 
saturation. Hence even with well resolved observations 
of the Na D lines it is impossible to derive a precise red- 
dening from these lines and published estimates should 
be viewed with caution. 

We have examined the implied outburst and quiescent 
luminosities, and the nature of the companion, in the 
context of these proposed scenarios. Neither a distance 
of 6 kpc, nor 15 kpc, presents difficulties in these respects. 
In fact the larger distance simplifies matters somewhat, 
for example in explaining the non-detection of the com- 
panion (likely a late-type sub-giant) in quiescence. Con- 
sequently we propose that a large distance, ~ 15 kpc, 
should be given consideration. 



We are grateful to Jeroen Homan, Jon Miller, and Seth 
Redfield for productive discussions on this topic and to 
Tom Maccarone for helpful comments on the manuscript. 
RIH is supported by NASA through Hubble Fellowship 
grant #HF-01150.01-A awarded by STScI. DS acknowl- 
edges a Smithsonian Astrophysical Observatory Clay Fel- 
lowship. JC acknowledges support from the Spanish 
Ministry of Science and Technology through the project 
AYA2002-03570. This work uses observations collected 
at ESO in Chile and the Anglo- Australian Telescope. We 
would particularly like to thank the ESO Director's Of- 
fice and VLT staff for a generous and efficiently executed 
award of Director's Discretionary Time. This work has 
also made use of the NASA ADS Abstract Service and we 
are grateful to Craig Markwardt for making the mpfit 
IDL fitting routines publicly available. 



Gallo, E., Corbel, S., Fender, R. P., Maccarone, T. J., & Tzioumis, 

A. K. 2004, MNRAS, 347, L52 
Gallo, E., Fender, R., k. Corbel, S. 2003, ATEL #196 
Garcia, M. R., McClintock, J. E., Narayan, R., Callanan, P., 

Barret, D., & Murray, S. S. 2001, ApJ, 553, L47 
Gelino, D. M., Harrison, T. E., & Orosz, J. A. 2001, AJ, 122, 2668 
Gray, D. F. 1992, The Observation and Analysis of Stellar 

Photospheres, Cambridge University Press 
Hameury, J.-M., Barret, D., Lasota, J. -P., McClintock, J. E., 

Menou, K., Motch, C, Olive, J.-F., & Webb, N. 2003, A&A, 

399, 631 
Home, K. 1986, PASP, 98, 609 

Hynes, R. I., Steeghs, D., Casares, J., Charles, P. A., & O'Brien, 
K. 2003, ApJ, 583, L95 



8 



Hynes et al. 



Fryer, C. L. & Kalogera, V. 2001, ApJ, 554, 548 

Kolb, U., King, A. R., & Baraffe, I. 2001, MNRAS, 321, 544 

Kong, A. K. H., McClintock, J. E., Garcia, M. R., Murray, S. S., 

Barret, D. 2002, ApJ, 570, 277 
Kuulkcrs, E., den Hartog, P. R., in't Zand, J. J. M., Vcrbunt, 

F. W. M., Harris, W. E., & Cocchi, M. 2003, A&A, 399, 663 
Maccarone, T. J. 2003, A&A, 409, 697 

McClintock, J. E., Narayan, R., Garcia, M. R., Orosz, J. A., 

Remillard, R. A., & Murray, S. S. 2003, ApJ, 593, 435 
Mauche, C. W. & Gorenstein, P. 1986, ApJ, 302, 371 
Miller, J. M. ct al. 2004, ApJ, 601, 450 
Munari, U. & Zwitter, T. 1997, A&A, 318, 269 
Nakanishi, H. & Sofue, Y. 2003, PASJ, 55, 191 
Orosz, J. A. & Bailyn, C. D. 1997, ApJ, 477, 876 
Predehl, P., Braeuninger, H., Burkert, W., & Schmitt, J. H. M. M. 

1991, A&A, 246, L40 
Press, W. H., Teukolsky, S. A., Vcttcrlink, W. T., Flannery, B. P., 

1988, Numerical Recipes in C, 2nd Edn., Cambridge University 

Press 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Sembach, K. R., Danks, A. C, & Savage, B. D. 1993, A&AS, 100, 
107 



Sembach, K. R. & Danks, A. C. 1994, A&A, 289, 539 
Shahbaz, T. 2003, MNRAS, 339, 1031 

Shahbaz, T., Fender, R., & Charles, P. A. 2001, A&A, 376, L17 
Shahbaz, T., van dcr Hooft, F., Casares, J., Charles, P. A., & van 

Paradijs, J. 1999, MNRAS, 306, 89 
Soria, R., Wickramasinghc, D. T., Hunstead, R. W., & Wu, K. 

1998, ApJ, 495, L95 
Soria, R., Wu, K., & Johnston, H. M. 1999, MNRAS, 310, 71 
Spitzcr, L. 1978, Physical Processes in the Interstellar Medium, 

Wiley, New York 
Sutaria, F. K. ct al. 2002, A&A, 391, 993 
Tomsick, J. A. ct al. 2003, ApJ, 597, L133 

Wagner, R. M., Starrfield, S. G., Hjellming, R. M., Howell, S. B., 

Kreidl, T. J. 1994, ApJ, 429, L25 
Welty, D. E. & Hobbs, L. M. 2001, ApJS, 133, 345 
White, N. E. & van Paradijs, J. 1996, ApJ, 473, L25 
Wu, K., Soria, R., Hunstead, R. W., & Johnston, H. M. 2001, 

MNRAS, 320, 177 
Zdziarski, A. A., Poutanen, J., Mikolajewska, J., Gierlinski, M., 

Ebisawa, K., & Johnson, W. N. 1998, MNRAS, 301, 435 



